Real-space fluctuations of effective exchange integrals in high-T^ 
cuprates 
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Abstract. - We present ab-initio calculations of effective magnetic exchange, J, as well as Hub- 
bard parameters {t, U and Se) as a function of the local distribution of doping atoms for the 
high- Tc superconducting (CaxLai_x)(Bai.75_xLao.25+x)Cu30y family. We found that both the 
exchange and the energies of the magnetic orbitals are strongly dependant on the local dopant 
distribution, both through the induced modification of the apical oxygen location and of the in- 
duced local electrostatic potential. The J real-space map, for a random distribution of dopants, 
positively compares with observed STS gap inhomogeneity maps. Similarly, the orbital energy 
fluctuations induce weak charge inhomogeneities on copper sites, that can be positively compared 
with the observed LDOS inhomogeneities. These results clearly support an extrinsic origin of both 
the gap inhomogeneities and LDOS. 
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Since their discovery by Miiller in 1986 [1], high-Tc su- 
perconductors have been the subject of an extensive re- 
search activity, both experimentally and theoretically. Re- 
cently, the question of the spatial (in)homogeneities of 
the superconducting state has attracted a lot of atten- 
tion. Indeed, in the recent years, the improvement and 
development of the tunneling microscopy techniques such 
as Scanning Tunneling Spectroscopy (STS) have allowed 
the access to spatial imaging of quantities such as the den- 
sity of states and the superconducting gap. Using these 
techniques, nanoscale spatial inhomogeneities of the gap 
amplitude were observed [2]. Strong fluctuations of the 
superconducting gap or even suppression of the supercon- 
ductivity was found within islands of typical size of one 
nanometer in the Bi2Sr2CaCu208+5 (Bi-2212) compound 
at 4.2K, that is far below the critical temperature [3]. The 
origin of these gap inhomogeneities is quite controversial, 
however crucial since some models link the high value of 
Tc to electronic phases separation or disorder. The con- 
troversy opposes authors defending the idea that this is 
an intrinsic feature of the superconducting state [3-5] , to 
authors defending the idea of an extrinsic origin [6] , linked 
to the chemical inhomogeneities. The first issue is directly 
related to the models originating the high- Tc superconduc- 
tivity in electronic phases separation [7]. Indeed, in such a 
case, the gap inhomogeneities would be a pure electronic 



effect in an otherwise (chemically and structurally) ho- 
mogeneous system. The second issue emphasizes the role 
of the always present chemical disorder resulting from the 
doping process out of the Mott-Hubbard insulating parent 
system. 

Despite the unknown origin of high- Tc superconductiv- 
ity it is widely accepted that the antiferromagnetic cor- 
relations play an important role in its onset, either as a 
competing phase [8] or as a pairing mechanism [9]. Some 
authors even claim a direct scaling of the critical temper- 
ature with the effective magnetic exchange, J, between 
adjacent copper atoms [10, 11]. In any case J plays a 
central role in the superconductivity onset and thus its 
fluctuations as a function of the chemical and structural 
inhomogeneities is expected to yield valuable new insight 
on the origin (intrinsic or extrinsic) of the observed gap 
and Tc spatial fluctuations. 

The present work proposes to study the local fluctu- 
ations of the effective exchange integral, J, as a func- 
tion of the chemical and structural disorder in the 
(CaxLai_x)(Bai.75_xLao.25-i-x)Cu30y (CLBLCO) copper 
oxides family. Indeed, this family presents two determi- 
nant characteristics, (i) the effect of hole doping and of 
chemical disorder is decoupled (ii) structural informa- 



simple formal charges analysis shows that the doping in the 
CLBLCO family is independent of the amount of calcium, x, in- 
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tions on the local distortions as a function of the chemical 
disorder are available [12]. 

For this purpose we used the CAS [13]+DDCI [14] 
(Complete Active Space + Difference Dedicated Con- 
figurations Interaction) quantum chemical spectroscopy 
method on embedded fragments. Indeed, this method is 
presently the most reliable and accurate one for the de- 
termination of effective magnetic integrals. It proved its 
efficiency on copper oxides superconductors and predicted 
the J values within experimental accuracy for a large num- 
ber of such compounds [10]. 

The CAS+DDCI method is an exact diagonalization 
method within a selected set of electronic configurations, 
specifically chosen so that to properly treat (i) the strong 
electronic correlation within the set of magnetic orbitals 
and (ii) the screening effects on the all magnetic configu- 
rations. The orbital space is thus divided into three sub- 
sets : the occupied orbitals set spanning orbitals al- 
ways doubly-occupied in the reference configurations (the 
CAS), the active orbitals set for which all possible or- 
bital occupations and spins are found in the CAS config- 
urations and the virtual orbitals set spanning orbitals 
always remaining empty in the CAS configurations. If 
the active orbitals are defined as the magnetic orbitals of 
the copper atoms, the CAS definition will thus insure the 
proper treatment (by exact diagonalization) of the corre- 
lation effects within them. The screening effects are then 
treated by adding to the diagonalization space all single 
and double excitations, on all CAS configurations, that 
participate to the excitations energies at the second order 
of the quasi-degenerate perturbation theory. 

The embedded fragments are divided into a quantum 
part and a bath. The quantum part includes the magnetic 
atoms (copper atoms), the ligands mediating the interac- 
tion as well as the first coordination shell of magnetic and 
ligands fragments. The bath is built in order to reproduce 
the main effects of the rest of the crystal on the quantum 
part, that is the Madelung potential and the exclusion ef- 
fects. The latter corresponds to the orthogonalisation of 
the quantum fragment orbitals to the orbitals of the rest 
of the crystal. This effect is usually taken into account 
through total ions pseudo potentials [15] (TIPs) located 
at the crystallographic positions of the first coordination 
shell of the quantum fragment. The former is reproduced 
using a set of renormalized charges. 

The (CaxLai„x)(I3ai.75„xLao.25+x)Cu30y compounds 
with 0.1 < a; < 0.4 and 6.8 < ?; < 7.25 are high tem- 
perature double-layers superconductors belonging to the 
YBCO group. As in the layered YBCO structure, the 
CLBLCO compounds are built from two different types 



troduced in the compound ; the hole doping only depends on the 
oxygen concentration y. Two different copper oxide planes are crys- 
tallographically present in CLBLCO family : the superconducting 
ones, associated with crystallographic Cu II sites, that do not see 
much fluctuations of their oxygen composition and the non super- 
conducting ones, associated with Cu I sites, with a, y — 6 oxygen 
composition. 



of copper planes : (i) two Cu02 superconducting lay- 
ers corresponding to the Cull crystallographic sites (see 
ref. [12]) separated by randomly distributed calcium and 
lanthanum ions (named Lay), and (ii) basal CuOy_6 lay- 
ers, corresponding to the crystallographic Cu I coppers. In 
the latter planes, the oxygen atoms (named Oq,) are ran- 
domly distributed along the Cu-Cu bonds. Let us note 
that on the contrary to the YBCO compounds, no ex- 
tended CuO chains are observed in these Cul layers. The 
YBCO Ba sites are randomly occupied with barium and 
lanthanum ions (named LaBa)- Associated with the apical 
oxygen atoms (Oc), they form layers in between the Cull 
and Cul planes. 

From the above crystallographic data [12] one can spot 
four origins for the chemical and structural disorder, (i) 
The random distribution of the Ba and Lasa ions in the 
YBCO Ba site, (ii) The in plane displacements of the Oc 
(apical) and Oq (Cul layers) oxygen atoms as a function 
of the doping ion on the YBCO Ba site. Indeed, the Ba'^^ 
ions are 20% larger than the La'^+ ones, (iii) The random 
distribution of Ca and Lay ions in the YBCO Y site, (iv) 
The random distribution of the Oa ions. 

It is usually accepted that the exchange couplings are 
mainly sensitive to the geometry of the coppers first coor- 
dination shells. In the CLBLCO family this geometry is 
controlled by the YBCO Ba sites occupation [12]. In the 
present work we thus focused on the effect of this degree 
of freedom, and averaged the electrostatic effects induced 
by the Oq, positional disorder as well as those induced by 
the Ca / Lay chemical disorder. 

We computed the effective exchange J between nearest 
neighbor (NN) copper atoms as well as the parameters of 
a Hubbard model (that is t, U and fc) as a function of the 
local ion distribution on the YBCO Ba site. Following the 
above specifications for building the embedded fragment, 
the quantum part was reduced to a CU2O9 fragment sur- 
rounded by 6 Cu2+ and 10 O^- tIPs in the Cu02 plane, 
2 Cul TIPs on top of the apical oxygens, 6 (Ca/La)(3-'^)+ 
averaged TIPs on the YBCO Y sites, and 6 either Ba^+ or 
La'^^ on the YBCO Ba sites, as pictured in figure [TJ Dis- 
tributing Ba^+ and La'^"'" ions in all possible manners in 
the 6 available positions next to the CU2O9 fragment, one 
finds 24 non-equivalent configurations. The probability of 
each distribution depends on the value of x. According to 
the X-Ray structural data of reference [12], each of these 
configurations is associated with a specific displacement 
of the apical oxygens of the CU2O9 fragment as well as a 
specific localization of the Ba^+ and La'^+ ions, accord- 
ing to their nature. For each configuration, the atomic 
positions used in our calculations (both for the CU2O9 
fragment and the TIPs) were extracted from the exper- 
imental data of reference [12]. The set of renormalized 
charges was built so that to reproduce the "long" range 
part of the Madelung potential of the crystal, averaged 
over ions disorder. The renormalization factors from the 
formal charges were computed along the method of refer- 
ence [16]. We imposed the nullity of all multipolar mo- 
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Fig. 1: (Color online) Schematic representation of the quantum 
part and TIPs (dashed) of the embedded fragments. The dif- 
ference between the quantum part and the TIPs is emphasized 
by used of undashed/dashed atoms and large/small bond lines. 

ments up to order 8. It resulted a set of 90616 renormal- 
ized charges. The charges locations were set to the exper- 
imental crystallographic positions given in reference [12], 
without the local displacements associated with a specific 
local ions configuration. The ions nominal charges were 
taken as the averaged nominal charges on the crystallo- 
graphic site ; that is for the superconducting Cu02 layers 
and related apical oxygens : Cu = 2+, O = 2—, for the 
doping CuOy_6 layers : Cu = 2-f , O = x 2— , where r] is 
the site occupation probability in the AFM parent system, 
for the Ba/La and Ca/La layers 2.125 -I- x/2 and 3 — x. 
The Madelung potential obtained from this set of renor- 
malized charges exhibited an error smaller than 1 mcv on 
any site of the quantum part. We also ran test calcula- 
tions in order to verify that (i) the charge modifications 
associated with the doping from the AFM parent, (ii) as 
well as the actual ion distribution in those far away layers, 
do not affect the results in any significant way (error bar 
2-3 meV). 

We computed the J, t, U and Se parameters for the 
24 distributions, for x = 0.1 and x — 0.4 and the hole 
doping parameter y associated with the maximal critical 
temperature given in reference [12] — that is y = 7.158, 
= 52.6 K for x = 0.1 and y = 7.174, = 80.3 K for 
X = 0.4. 

Figure [2] presents the probability of each J value occur- 
rence and figure [3] displays the resulting real-space map of 
magnetic exchange obtained for a random distribution of 
the Ba and La ions in the YBCO Ba sites. One sees im- 
mediately that the J values are spread over a range wider 
than 20 meV for each system despite the fact the Cu- 
0-Cu buckling angle does not depend on the Ba / Lasa 
distribution. Indeed, the buckling angle is only system 
dependant. For a given compound, the J distribution is 
thus independent of it, while this is the contrary for the 
average exchange value, J. One finds, for the two studied 
compounds, Jx=o.i = — llOmeV for a buckling angle of 
e = 167.2° and J^=o.4 = -132 meV for = 170.0°. Let 
us note that the variation of J with 9 is as expected, that 
is the super-exchange term increases with the overlap be- 
tween the magnetic Cu 3dx'^_y2 orbitals and the bridging 
in-plane oxygen 2p orbitals (i.e. with cos^ {6 — tt)). 
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Fig. 2: (Color online) Probability of effective exchange val- 
ues. Antiferromagnetic exchange J are taken negative. The 
inset shows the probability of occurrence for effective ex- 
change values when the atomic displacements (induced by 
the nature of the surrounding Ba site ions) of the apical 
oxygens are ignored. Solid blue lines are associated with 
(Cao.iLao.9)(Bai.65Lao.35)Cu307.i58 and dashed red lines with 
(Cao.4Lao.6)(Bai.35Lao.65)Cu307.i74. 
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Fig. 3: (Color online) Real-space map of magnetic exchange 
obtained for a random distribution of the Ba and La ions in 
the YBCO Ba sites in the (Cao.4Lao.6)(Bai.35Lao.65)Cu307.i74 
compound. 



Let us now analyze the possible origin of the fluctuations 
around these average values. The buckling angle and the 
embedding can be eliminated since they are unique for 
each compound. The only remaining degrees of freedom 
are related with the local Ba / LaBa distributions. One can 
think of i) the modification of the apical oxygen (Oc) lo- 
calization as a function of the local Ba / Lasa distribution, 
and ii) the electrostatic potential fluctuations induced by 
the different Ba^"*" / La^"^ charges. The relative influence 
of these two degrees of freedom can easily be determined 
by computing the J values for each Ba / Lasa distributions 
but with undistorted geometries (see inset of figure^]). Do- 
ing so, one sees that i) the J values still present a large 
dispersion and ii) that the dispersion range is much larger 
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(nearly double) than when the apical oxygen displace- 
ments arc taken into account. It results that contrarily 
to what is usually admitted, the electrostatic effects play 
an important role in the effective exchange values. In fact, 
one of the mayor role of the atomic position relaxation is 
to moderate this effect. One can thus expect that the ac- 
tual spatial fluctuations of J is even larger than the ones 
calculated in this work and displayed in figure [3l Indeed, 
not only the electrostatic disorder induced by the Ca / La 
distribution on the YBCO Y site has not been explicitely 
taken into account in the present calculations, but in ad- 
dition there are not any atomic displacement associated 
with these Ca / La distributions that could soften their 
electrostatic effects. 
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Fig. 4: (Color online) Probability of effective transfer 
integrals values. Solid blue lines are associated with 
(Cao.iLao.9)(Bai.65Lao.35)Cu307.i58 and dashed red lines with 
(Cao.4Lao.6)(Bai.35Lao.65)Cu307.i74. 

Figure 2] displays the probability of occurrence of the 
effective hopping values between NN copper atoms. The 
average transfer for the two compounds t is found to be 
tx=o.i = — 493meV and tx=oA = — 537meV, while the 
ranges of variation are respectively 36meV and 37meV. 
Basically, the effective transfer integrals behave in a simi- 
lar way as the effective exchange, both as a function of the 
buckling angle and of the local electrostatic fluctuations. 
As far as the on-site repulsion U is concerned, the sys- 
tem average values are very similar with Ux=q.i — 8.9 eV 
and Ux=QA = 8.8 eV, and the spatial fluctuations very 
small in relative values with {AU/U)x=o.i — 0.07 and 
{AU/U)x=o.4 = 0.06. Let us note that the J fluctuations 
nicely agree with their perturbative evaluations from t and 
U fluctuations, that is A J ~ -8{t/U)At + A{t/U)'^AU. 

From the ab initio calculations (excitation spectrum and 
associated wave functions) one can also derive the differ- 
ence between nearest neighbor magnetic orbital energies : 
ds. As for the exchange, transfer and on-site Coulomb re- 
pulsion the effective energies of the magnetic orbitals do 
strongly fluctuate in real space. In fact, the orbital en- 
ergy difference between NN sites is proportional to the 
difference in the number of Ba / La ions surrounding each 
copper atom (see figure [5]). It comes an energy stabiliza- 
tion of 1.4 eV per additional Ba La substitution. One 
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Fig. 5: Magnetic orbital energy difference between NN sites as 
a function of the charge difference associated with the coppers 
first neighbors YBCO Ba sites. Crosses: computed points, 
dashed line: linear fit. 
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Fig. 6: (Color online) Real-space map of the number of Ba / La 
ions surrounding a given copper site (or equivalently of the 
magnetic orbital energy) for a random distribution of the 
Ba / La atoms in the YBCO Ba sites. 



can thus expect an associated real-space electronic density 
fluctuation as a function of the local ionic conflguration. 
Indeed, our calculations show a correlated fluctuation in 
the copper atoms population of about 0.06e. Figure [B] 
pictures such an expected fluctuation for a random distri- 
bution of the Ba / La atoms in the YBCO Ba sites. The 
typical scale of the magnetic orbital energies fluctuations 
and thus of the expected density of states is of the order 
of a few nanometers. 

In summary our calculations showed that the extrinsic 
inhomogeneities associated with material doping is respon- 
sible for both magnetic exchange and on-site orbital ener- 
gies strong spatial fluctuations. In both cases, the fluctu- 
ations are of a typical scale of a few nanometers. First, 
the computed magnetic orbital energy fluctuations induce 
an associated fluctuation of the copper density of states 
and positively correlate (similar typical size and pattern) 
with the observed LDOS maps as picture in figure lb of 
reference [3]. Second, the strong relation between mag- 
netism and superconductivity in high-Tc cuprates [8-11] 
induces us to compare our exchange fluctuations with the 
gap inhomogeneities observed in STS experiments [2,3]. 
Despite the fact that it relates to different compounds, the 
typical size of doping inhomogeneities is similar and flg- 
ure[2]of the present work compares positively with, for in- 
stance, figure 1 of reference [4], in all aspects except for the 
amplitude of the inhomogeneities. However as discussed 
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above, the electrostatic disorder of the YCBO Y sites was 
not taken into account and slioukl strongly increases this 
amplitude. As a matter of example, we computed the J 
modification associated with a strongly dissymmetric dis- 
tribution of the Ca / La ions on the YCBO Y sites and 
found that the exchange integral is being changed from 
-146 meV to -98meV. 

In the Bi2Sr2CaCu208+<5 compounds, the chemical and 
structural inhomogeneities located close to the supercon- 
ducting layers are, of course, not due to the counter-ions 
(as in the present compounds) but rather to the non sto- 
ichiometric oxygen atoms. Our results clearly show that 
this is i) the fluctuations of the apical oxygen location 
and ii) the local electrostatic environment of the supercon- 
ducting copper plane that are responsible of the exchange 
and orbital energy fluctuations. The chemical origin of 
these structural and electrostatic effects are system depen- 
dent, related to doping counter-ions in our case, related to 
the non stoichiometric oxygen in Bi2212 compounds. In 
any case our results support the conclusions of McElroy et 
al [17] of an extrinsic (chemical and structural) origin of 
the gap inhomogeneities and associated critical tempera- 
ture. 
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